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Application  of  the  Ring  Vortex  Method  to  Aerodynamic  Design  of  Lifting  Rotor 
Systems 

By:  V,  I.  Shaydakov 
Introduction: 

The  operation  of  the  lifting  rotor  during  flight  with  a  horizontal  velocity 
component  has  been  the  object  of  continuous  studies  by  a  number  of  foreign  and 
Soviet  scientists,  particularly  in  the  past  decade*  At  the  present  time,  great 
successes  have  been  achieved  in  the  area  of  determination  of  the  forces  and 
moments  acting  on  the  lifting  rotor  in  a  slanted  stream,  on  the  basis  of  both 
flat  and  three-dimensional  models  of  the  vortex  shroud. 

However,  the  conclusions  of  the  existing  theory  cannot  be  extended  to  the 
mode  of  rapid  descent  of  a  helicopter  with  and  without  application  of  power  to 
the  lifting  rotor  shaft.  These  are  the  vortex  ring  and  turbulent  (autorotation) 
modes.  As  we  know,  the  stream  theory  of  Glauert  [1]  produces  results  for  these 
modes  which  do  not  agree  with  experimental  results.  One  frequent  case  of  vortex 
ring  and  autorotation  modes  is  vertical  descent  of  a  helicopter,  which  has  been 
studied  by  the  author  of  the  present  article  on  the  basis  of  general  theories 
of  aeromechanics  [5]  and  by  Ye.  S.  Vozhdayev  on  the  basis  of  vortex  theory. 

Steep  glide  modes  of  a  helicopter  correspond  to  an  extremely  complex  pic¬ 
ture  of  flow  around  the  lifting  rotor,  which  makes  experimentation  difficult  and 
does  not  allow  the  vortex  model  of  the  phenomenon  to  be  constructed. 

All  of  this  led  us  to  the  idea  of  simplifying  the  existing  vortex  models 
of  the  lifting  rotor  as  much  as  possible  and  finding  the  simplest  method  of 
determination  of  its  aerodynamic  characteristics,  in  order  to  use  this  method 
later  to  investigate  the  modes  of  steep  descent  of  a  helicopter  in  the  first 
approximation . 

This  article  presents  a  method  (which  we  will  refer  to  as  the  ring  vortex 
method) ,  providing  the  same  accuracy  as  the  theory  of  the  ideal  rotor  or  the 
vortex  NEZH  rotor  theory. 

1,  Content  of  suggestion  method 

N.  Ye.  Zhukovskiy  proved  a  theorem  [4]  for  determination  of  the  resistance 
of  a  flat  circular  plate  placed  perpendicular  to  a  flow,  which  can  be  formulated 
as  follows:  the  momentum  (M)  imparted  by  a  closed,  flat  vortex  filament  to  a 
limitless  liquid  mass  in  the  direction  perpendicular  to  the  plane  in  which  it 
is  located  is  equal  to  the  product  of  the  density  of  the  fluid  (p)  times  the 
circulation  of  the  vortex  (r)  and  the  area  which  limits  it  (F)  (Figure  1) 

(1] 
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It  has  been  determined  that  this  theorem  can  be  used  in  the  solution  of 
certain  problems  of  aerodynamics.  If  we  know  how  many  closed  vortices  are 
formed  per  unit  time  or  the  extent  to  which  the  area  F  limited  by  them  in- 
crease'?  ne.-r  unit  timf^  can  calculate  the  change  in  the  momentum  and 
conse(  :ting  on  the  system. 


Figure  1, 


Figure  2. 


To  illustrate  this,  let  us  analyze  the  case  of  a  wing  of  finite  length 
with  a  constant  circulation  along  the  span.  Figure  2  shows  the  vortex  system 
of  a  wing  moving  at  velocity  V^.  The  attached,  acceleration  and  free  vortices 

make  one  closed  vortex,  the  plane  of  which  is  inclined  to  the  direction  of  the 
unperturbed  flow  by  angle  e  [e  is  the  inductive  taper  angle  of  the  flow] .  The 
lift  R  is  perpendicular  to  the  plane  of  the  closed  vortex  and  produces  a  com¬ 
ponent  in  the  direction  of  the  velocity  of  the  unperturbed  flow  (inductive 
drag).  We  will  base  our  considerations  on  an  ideal  fluid,  in  which  case  this 
vortex  system  can  exist  for  an  infinite  time. 


The  motion  of  the  wing  is  accompanied  by  an  increase  in  area  F,  limited 
by  the  closed  vortex,  causing  a  change  in  the  momentum.  For  example,  over 
time  At  the  wing  travels  a  distance  of 


Kp 


o» 


and  the  free  vortices  are  elongated  by  As^^: 

Ld 


As, 


-An',. 


cos.  cost 

In  this  same  time,  a  change  in  the  momentum  by  AM  occurs: 

AM  =  pFAF-  pFA^^.-/  = 
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The  force  R  in  one  second  is  equal  to  the  change  in  momentum  per  second: 

/?  =  ^  =  pri^,/.  (2) 

This  formula  is  known  as  the  Zhukovskiy  formula  for  the  lift  of  a  wing. 

It  is  particularly  convenient  to  use  the  method  suggested  to  determine 
the  aerodynamic  characteristics  of  the  lifting  rotor  of  a  helicopter  in  axial 
and  inclined  streams, 

2.  Statement  of  Problem.  Thrust  Problem 

Let  us  study  a  NEZH  rotor  with  constant  circulation  r  over  the  radius. 

A  free  vortex  leaves  the  tip  of  each  blade  with  the  same  circulation  P 
enclosed  in  the  ring  as  soon  as  the  blade  completes  one  rotation.  In  one 
rotation,  a  rotor  with  k  blades  emits  k  ring  vortices,  limiting  area  ttR^  , 
equal  to  the  area  of  the  rotor  disk. 

Thus,  an  operating  rotor  forms  a  vortex  track  in  its  wake,  consisting  of 
discrete  ring  vortices,  moving  in  the  direction  of  the  stream  striking  the 
rotor.  The  vortices  have  identical  circulation  F  and  are  located  in  planes 
parallel  to  the  plane  of  rotation  of  the  lifting  rotor.  In  an  ideal  fluid  the 
vortices  do  not  break  up;  therefore,  the  vortex  track  has  infinite  length. 

In  the  general  case,  the  ring  vortices,  moving  over  the  flow  surface,  change 
their  dimensions  due  to  the  presence  of  constriction  of  the  stream  following 
the  rotor. 

Figure  3  shows  the  vortex  model  which  we  have  described  for  the  case  of 
operation  of  a  rotor  in  an  inclined  flow. 


Let  us  now  go  over  to  calculation  of  the  thrust  of  the  lifting  rotor  (T) . 
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The  time  of  formation  of  one  ring  will  be 

CO  is  the  angular  velocity  of  rotation  of  the  rotor. 


In  time  At,  momentum  AM  is  created,  equal  to  the  force  T: 

Considering  (3) ,  we  have:  ■ 


(3) 


(4) 


As  we  can  see,  (4)  is  suitable  for  any  mode  of  operation  of  the  rotor.  This 
same  formula,  although  somewhat  longer,  can  be  produced  fi'om  the  theory  of  an 
isolated  blade  element. 

We  must  now  find  the  relationship  between  circulatioii  T  and  the  inductive 
velocity  in  the  plane  of  the  rotor  disk.  To  do  this,  let  us  go  over  to  study 
of  the  vortex  track  behind  the  blade. 

The  formulas  produced  ar.e  known  and  indicate  that  the  method  suggested  is 
sound  in  its  basis,  simple  and  contains  no  grossly  incorrect  assumptions. 

In  the  next  section,  we  will  study  the  modes  of  a  lifting  rotor,  the 
mathematical  description  of  which  has  always  involved  significant  difficulty. 

3.  Rotor  in  Slanted  Flow  Mode 

The  first  theoretical  studies  of  these  modes  were  performed  in  1928  by 
Glauert,  who  studied  the  average  inductive  velocities  v^,  over  the  disk.  He 

based  his  calculations  on  the  hypothesis  that  the  lifting  rotor  in  operation 
acts  on  a  stream  of  air,  the  cross  section  of  which  has  a  diameter  equal  to 
the  diameter  of  the  lifting  rotor.  This  hypothesis  is  similar  to  the  assump¬ 
tion  of  Prandtl  for  a  wing  [1] . 

The  thrust  formula  in  this  case  will  be 


r  =  (5) 

where  is  the  rate  of  flow  of  air  through  the  lifting  rotor  disk.  This  for¬ 
mula  is  broadly  used  in  aerodynamic  calculations  involving  lifting  rotors  acting 
in  slanted  flows  [2,6]  and  for  determination  of  the  slope  of  the  vortex 
solenoid  (e.g.,  see  [3]). 
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The  vortex  ring  method  allows  us  to  give  a  simple  proof  for  Glauert's 
formula.  Figure  4  shows  the  vortex  track  resulting  from  operation  of  a 
rotor  in  the  slanted  flow  mode.  Let  us  separate  two  characteristic  cross 
sections.  Section  1-1  lies  in  the  plane  of  the  rotor  disk.  Here  the  air 
flow  rate  through  the  disk  is  represented  as  the  vector  sum 

K-Vo+it,  (6) 

where  v^  is  the  mean  inductive  velocity  in  cross  section  1-1,  resulting  from 
the  vortex  track. 

Section  2-2  lies  far  behind  the  rotor.  Here  the  stream  has  undergone 
maximum  constriction,  and  the  vortex  track  has  been  converted  to  an  inclined 
cylindrical  vortex  solenoid.  The  velocity  within  the  solenoid 

^2=^0  +  ®2,  (7) 


where 


is  the  mean  inductive  velocit/  resulting  from  the  vortex  solenoid. 


Figure  4. 

It  can  be  considered  with  a  high  degree  of  accuracy  that  the  vortex  rings 
follow  each  other  closely,  forming  a  continuous  vortex  layer.  Then  the  in¬ 
ductive  velocity  within  the  vortex  track  increases  continually  and  monotoni- 
cally  with  increasing  distance  from  the  rotor,  reaching  its  maximum  at  cross 
section  2-2. 

Outside  the  vortex  solenoid  in  cross  section  2-2,  the  mean  inductive  ve¬ 
locity  is  equal  to  zero,  and  the  total  velocity  is  equal  to  the  velocity  of 
the  unperturbed  stream  V^.  The  mean  running  circulation  of  the  vortex  layer 

around  the  circumference  in  the  direction  of  the  axis  of  the  inclined  vortex 
cylinder 
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(8) 


Tb.c  ’^1^2-  i^oI  =  l®2l, 


while  the  velocity  of  motion  of  the  vortex  layer  (i,e.,  the  velocity  of  motion 
of  the  ring  vortices)  is  the  vector  half  sum  between  the  velocities  and 

1  r  .  ^ 


Va  +  V^o. 
2 


(9) 


It  follows  from  the  Froud-Finsterwalder  theorem  [7]  that  for  all  operating 
modes  the  rotor 


1)2  =  2®!. 

This  theorem  is  called  the  theorem  of  velocities  in  the  stream  of  a  rotor  and 
is  proven  on  the  assumption  that  the  vortex  column  behind  the  rotor  is  not 
constricted.  (In  our  case^  the  vortex  rings  do  not  change  their  dimensions 
as  they  move  away  from  the  rotor.)  In  the  practice  of  aerodynamic  calculation, 
constriction  of  the  stream  beyond  the  rotor  is  generally  ignored.  This  pro¬ 
duces  no  great  divergence  with  experimental  results.  Therefore,  we  can  consider 
that 


and  _  ^  ^ 

i>..c  =  Vo  +  vi^y,. 


(11) 

(12) 


The  running  circulation  of  the  vortex  layer  in  Section  2-2  can  be  produced 
by  another  method  as  well . 

During  the  time  of  formation  of  one  ring  (At)  the  vortex  ring  in  Section  2-2 
moves  over  a  path  of 


As  =  A/.  V,.c  . 


(13) 


Substituting  formulas  (3)  and  (12)  here,  we  produce 


(14) 


The  running  circulation  of  the  vortex  layer  in  the  direction  of  the  axis  of  the 
vortex  track  will  be; 

_  r  _  *ra> 

1b.  C  ~As“2jtV',‘  .  (18) 


Setting  equations  (11)  and  (15)  equal  to  each  other,  we  find: 
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AFu)  =  4111^1111. 


C16) 


Substituting  (16)  into  formula  (4)  for  the  thrust,  we  find: 


T  = 


(17) 


Here  m  is  the  mass  flow  rate  of  air  through  the  disk  of  the  lifting  rotor  per 
second.  As  we  can  see,  the  formula  produced  corresponds  fully  with  (5),  which 
was  to  be  proven.  We  will  relate  all  velocities  to  the  inductive  velocity  of 
hovering  (v^^): 


■ 


Then  formula  (17)  becomes: 


1  = 


(18) 


It  follows  from  Figure  4  thafl 

'^i  Vo  +  Vi  —  2^,d, sln-r. 


(19) 


Let  us  substitute  expression  (19)  into  (ISJ: 


rJ-2VoSlnx^f  f  1  =0. 

The  equation  produced  gives  us  the  relationship  between  inductive  velocity, 
the  velocity  of  the  oncoming  stream  and  angle  of  attack  t.  In  [6]  we  find 

the  graph  of 

We  can  attempt  to  estimate  the  influence  of  the  nonlinearity  of  the  axis 
of  the  vortex  cylinder  on  the  inductive  characteristics  of  the  rotor. 

Actually,  it  must  be  considered  that  in  cross  section  2-2,  the  vortex 

rings  continue  their  motion  in  the  direction  of  velocity  not  ,  as  occurs 

in  cross  section  1-1.  In  this  case,  the  rate  of  displacement  of  the  vortices 

in  cross  section  2-2  will  be  (Figure  4): 


^In  this  and  all  subsequent  formulas,  angle  r  should  be  substituted  with  the 
proper  sign,  in  our  case  with  the  minus  sign.  (The  angles  and  velocities  are 
taken  according  to  B.  N.  Yuriev  [6].) 


Vi.  C  =  cos  e. 


(21) 


where  e  is  the  angle  between  axes  of  the  vortex  solenoid  in  cross  sections 
1-1  and  2-2  or  the  angle  between  velocities  and  ^2' 

Then  formula  (18)  becomes 


1  =  V/j  Uj  cos  e. 


(22) 


while  equation  (20)  is  rewritten  as: 

+  - L_  =  o. 

COS  c 


(23) 


As  we  can  see, 
ductive  velocities. 


consideration  of  nonlinearity  slightly  increases  the  in- 
Let  us  determine  how  significant  this  increase  is. 


It  follows  from  the  triangle  of  velocities  in  cross  section  2-2 
(Figure  4)  that 


cos  e  = 


or  after  simple  conversions 

cos  £  — -  -- 


Vl  -f  —  3 VoU,  sin  t 
-  2VQVi  srn  r)  4-  4v^  —  ‘  '^) 


(24) 


In  order  to  determine  must  solve  equations  (23)  and  (24) 

jointly.  The  solution  is  quite  cumbersome.  We  can  use  a  much  simpler  approach. 
As  will  be  demonstrated,  consideration  of  nonlinearity  of  the  axis  of  the  vortex 
track  changes  the  inductive  characteristics  of  the  lifting  rotor  very  little; 
therefore,  cos  c  will  be  calculated  using  formula  (24),  substituting  the  values 
of  =  f(V^,T)  preliminarily  found  from  equation  (20),  which  ignores  this  non¬ 
linearity,  into  formula  (24). 


The  greatest  influence  of  nonlinearity  of  the  stream  axis  beyond  the 
rotor  is  on  the  mode  of  a  flight  with  zero  angle  of  attach  (-t  =  0)  . 

We  will  estimate  this  influence  using  coefficient  which  shows  the 
number  of  times  by  which  the  inductive  velocity  calculated  according  to  equation 
(20)  must  be  increased  in  order  to  consider  this  nonlinearity. 

In  the  table  we  present  the  results  of  calculation  for  the  case  -t  =  0: 
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V, 

0 

0,25 

0,5 

0,75 

1.0 

1.25 

1.5 

1.75 

2,0 

2.5 

cos  e 

3.0 

0,980 

0,974 

j 

0.954 

0.942 

0,940 

0,950 

0,965 

0.977 

0.989 

uo 

1,003 

1,007 

1,014 

1.021 

1,024 

j 

1.022 

'1.017 

. 

1,012 

i 

1,006 

As  we  can  see  from  the  table^  the  maximum  influence  of  nonlinearity  of 
the  axis  of  the  vortex  cylinder  on  inductive  velocity  does  not  exceed  2,4%. 
Therefore,  the  nonlinearity  of  the  stream  axis  behind  the  rotor  can  be  ig¬ 
nored  in  vortex  theories,  which  is  usually  done. 

Thus,  the  method  we  have  presented  can  be  successfully  used  in  aero¬ 
dynamic  design  of  lifting  systems,  particularly  the  lifting  rotors  of  heli¬ 
copters.  The  replacement  of  the  vortex  cylinder  with  a  system  of  discrete 
vortex  rings  is  equivalent  to  expansion  of  the  vortex  cylinder  into  rings 
and  longitudinal  vortices,  which  is  always  performed  in  vortex  rotor  theories. 
Ignoring  the  longitudinal  vortices  is  equivalent  to  ignoring  the  closure  of 
the  stream  behind  the  rotor.  For  low-loaded  lifting  rotors,  as  wc  know,  no 
great  error  is  thus  produced. 

This  method  can  be  used  for  rotors  with  variable  circulation  over  their 
radius.  In  this  case,  integral  formulas  are  produced. 
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